Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 1 February 2008 



(MN MfeX style file vl.4) 



Dust and gas in luminous infrared galaxies - results from 
SCUBA observations 



U. Lisenfeld 1 , K.G. Isaak 2 * , R. Hills 3 

1 IRAM, Avenida Divina Pastora 7, N.C., 18012 Granada, Spain, e-mail: ute@iram.es 
2 Department of Astronomy, University of Maryland, College Park, USA 

3 Cavendish Astrophysics Group, Cavendish Laboratory, Madingley Road, Cambridge CBS OHE, UK 



OS 
<3\ 



1 February 2008 



o 
O 

00 
> 

in 
m 
O 

o 

0\ 



ABSTRACT 

We present new data taken at 850 /jm with SCUBA at the JCMT for a sample of 19 
luminous infrared galaxies. Fourteen galaxies were detected. We have used these data, 
together with fluxes at 25, 60 and 100 /im from IRAS, to model the dust emission. 
We find that the emission from most galaxies can be described by an optically thin, 
single temperature dust model with an exponent of the dust extinction coefficient 
(fc A oc A _/3 ) of /3 ~ 1.4 - 2. A lower ~ 1 is required to model the dust emission 
from two of the galaxies, Arp 220 and NGC 4418. We discuss various possibilities for 
this difference and conclude that the most likely is a high dust opacity. In addition, 
we compare the molecular gas mass derived from the dust emission, Msso^m, with the 
molecular gas mass derived from the CO emission, Mqo, and find that Mqo is on 
average a factor 2-3 higher than -Mgso^m- 
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1 INTRODUCTION 



Much observational evidence has been accumulated since 
the first discovery of luminous (infrared luminosities of 
Lir = L(8 - lOOO/Ltm) > 10 n L Q , "LIRGs") and ultralumi- 
nous (Lm > W 12 Lq, "ULIRGs") infrared galaxies to sug- 
gest that their interstellar medium (ISM) is quite different 
from that of normal galaxies. ULIRGs and LIRGs contain 
very large quantities of molecular g cLS, ELS traced by CO (e.g. 
Sanders et al. 1988), and by tracers of high density gas such 
as HCN and CS(3-2) (Solomon, Downes & Radford 1992). 
In spite of their abundance of molecular hydrogen, LIRGs 
have ratios of far-infrared (FIR) luminosity to molecular gas 
mass, Z/fir/ Mh, , that largely exce ed those measured in our 
own galaxy (Solomon et al. 1992). The ratio between the 
luminosity of the high-density tracers and Lfir, however, is 
quite normal. The FIR emission is a measure of the star for- 
mation activity in a galaxy and so the constant Lfir/Lhcn 
implies that the star formation rate (SFR) per dense gas 
mass is normal. The ratio of the luminosity of the high- 
density gas tracers to CO is much higher than in our own 
galaxy, and so it would appear that a larger frac tion o f the 
molecular gas is in a dense state. Solomon et al. (1997) con- 
clude from this that the ISM of ULIRGs is much denser than 
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in normal galaxies, with even gas in intercloud regions being 
in the molecular state. 

The dust properties of LIRGS are also different from 
those found in less luminous galaxies. There are a number 
of strong indications that the dust opacity in LIRGs is very 
high, with the result that the emission from LIRGs is opti- 
cally thick even in the FIR wavelength range. Based on the 
tight correlation between the flux at 100 um and the CO 
emission of ULIRGs, Solomon et al. (1997) argue that the 



ISM is optically thick at wavelengths as long as 100 fim. A 
high opacity in Arp 220 i s supp orted by a range of observa- 
tions: (i) Emerson et al. (1984) have concluded from obser- 
vations in the IR and submillimetre that the dust opacity 
is of order unity between 100 and 200 /im. (ii) From obser- 
vations with ISO usin g the Long Wavelength Spectrometer 
(LWS) Fischer et al. (1997) derived unity optical depth at 
~ 150/im. (iii) The ra tios of the IR fine structure line emis- 



sion observed by ISO (Genzel et al. 1998) are consistent with 
a screen opacity of Av — 45 mag, equivalent to an opacity 
of Av — 1000 mag if the dust and gas were well-mixed, 
(iv) Calculations by Downes & Solomon (1995) using mea- 
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surements of the CO surface brightness and a galactic dust- 
to-gas ratio suggest an opacity of about Ay ~ 1000 mag. 
The dust opacity is also high in NGC 4418, where a high 
dust extinction has been measured at optical wavelengths 
[Av 3> 50 mag, Roche et al. 1986). The high dust opacity 
might also explain why the C + 158 fim line is so weak in 
ULIRGs (Malhorta et al. 1997, Luhman et al. 1998). 
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An accurate measurement of the molecular gas mass 
of a galaxy is critical to determining star formation activ- 
ity. Estimating the molecular gas mass of an object is not 
straightforward, and is typically made in one of two ways: 
The first method relies on observations of the optically thick 
CO (1-0) emission. The calculated CO line luminosity can be 
converted to a H2 gas mass by adopting a conversion factor 
between CO and H2 that has been derived from observa- 
tions of Galactic molecu lar clouds (e.g. You ng & Scoville 
1982) and 7-rays fluxes (Bloemen et al. 1986j ). The second 
method makes use of the empirical relationship between gas 
and dust mass. The dust mass can be determined from the 
dust emission spectrum with assumptions made about the 
dust emissivity and composition. The dust mass can in turn 
be converted into a molecular gas mass using the locally 
derived gas-to-dust mass ratio (e.g. Hildebrand 1983). 

Both derivations assume that the dust and gas parame- 
ters derived from measurements in our own galaxy are valid 
under the quite different conditions in LIRGs, and that they 
are independent of star form ation activity. Theoretical argu- 
ments by Maloney & Black (1988) suggest that the CO-to- 
H2 conversion factor should be considerably lower in actively 
star forming galaxies because of the higher density and tem- 
perature of the constituent molecular clouds. Observation- 
ally, there are also indications that the Galactic CO-to-Hb 
conversion factor overestimate s the molecular gas mass in 
LIRGs: Shier, Rieke & Rieke (1994) find that the molecu- 
lar gas mass derived from CO emission using the standard 
Galactic conversion factor exceeds the dyna mical mass by a 
factor of 4 - 10. Similarly, Solomon et al. (1997) conclude 
from comparisons of the dynamical and gas mass derived 
from the FIR emission and the CO, that the Galactic con- 
version factor between the CO luminosity, Leo, and the 
molecular gas mass, Mh 2 , overestimates Mu 2 by a factor of 
about 5. This conclusion is supported by a better determi- 
nation of the dynamical mass of a sample of ULIRGs using 
high angular resolution observations made with the Plateau 



de Bure interferometer (Downes & Solomon 1998). 

In this paper we report on observations of the 850 
/im emission for a sample of 19 LIRGS galaxies using 
the submillimetre-wave continuum bolometer array, SCUBA 
( |Holland et al. 1999| ), at the JCMljj} The aim of these obser- 
vations has been twofold: firstly, to study the dust emission 
from a selection of galaxies believed to have an ISM quite 
different to that of our own Galaxy, and secondly, to com- 
pare the molecular gas mass derived from the dust emission 
at 850 /im with that derived from CO measurements in order 
to assess whether the assumptions on which both methods 
rely are valid for galaxies spanning a wide range of different 
star formation activities. The large field of view of SCUBA 
enables us to include any extended dust emission and to 
measure its size. In this way we are able to avoid the limita- 
tions of previous studies made with single-pixel bolometers 
(e.g. Andreani & Franceschini 1996, Lisenfeld et al. 1996). 



t The James Clerk Maxwell Telescope is operated by The Joint 
Astronomy Centre on behalf of the Particle Physics and Astron- 
omy Research Council of the United Kingdom, the Netherlands 
Organization for Scientific Research, and the National Research 
Council of Canada. 



2 SAMPLE AND OBSERVATIONS 

The sample consists of 19 galaxies selected from the sub- 
sample of the Bright IRAS galaxies observed in C O(l-O ) 
with the NRAO 12m telescope by Sanders et al. (1991). 
The combination of the large CO beam (55") and the small 
optical diameter of our sample means that no CO emis- 
sion has been missed. The galaxies were selected to cover 
a wide range of infrared luminosities (7.6 10 10 Lq < Lir < 

3.1 10 12 Lq) [j]and infrared-to-blue-luminosity ratios Lir/Lb 
(2 < Lir/Lb < 80), thus sampling a wide range of luminosi- 
ties and star formation activities (as traced by Lir/Lb). 

The galaxies were mapped in service time with SCUBA 
at the JCMT between July 1997 and April 1998. The 16- 
jiggle mode was used for the observations to obtain fully 
sampled maps at 850 ^m; the weather conditions were too 
poor in order to make use of the 450 /jm data that was ob- 
tained parallely. The chop throws ranged betweeen 60" and 
120" (see Table 1) and the integration time per map was 
between 15 and 70 minutes. The sky opacity at 850 /im var- 
ied between 0.2 and 0.45 . The data we re reduced using the 
software package SURF (Jenness 1998). The data were flat- 
fielded, extinction-corrected and despiked, with bad bolome- 
ters and bad channels removed prior to rebinning and flux 
calibration. Flux calibration was achieved using Uranus and 
the secondary calibrators CRL 618, IRC+10216, OH231.8. 
The conversion factor between volts and Jy was found to be 
quite consistent over different nights, with an average value 
of 264 ± 19 Jy/V. 



3 RESULTS 

Listed in Table 1 are the measured 850 (im fluxes (column 6) 
together with 12, 25, 60 and 100 jjm I RAS fluxes taken from 
the IRAS Faint Source Catalog ( jl990| ). The flux errors given 
in column 6 represent the noise in the maps, a n . The total 
error, er, used for the model fitting, takes additional account 
of the estimated calibration error, (j ca i, of an estimated 20 
per cent: a 2 = a 2 + a 2 ^. Fourteen of the nineteen galaxies 
were detected at 850 /jm, with 3a n upper limits of between 
50 and 200 mjy obtained for the remaining five galaxies. 
These upper limits have been computed under the assump- 
tion that the angular size of the galaxy is small compared to 
the beam; if considerable extended emission is present, the 
total fluxes could be higher. 

We have found extended emission for NGC 5653 and 
NGC 5936. For NGC 5653 a Gaussian fit yielded a Full 
Width Half Maximum (FWHM) of 17", whereas the shape 
of NGC 5936 is rather irregular with a maximum diameter 
of about 40". In the case of NGC 1614 and NGC 3110 there 
could be extended emission present, but the noise level of 
our data is not low enough in order to definitely confirm 
this. For the other galaxies we determined the source sizes 
to have a FWHM of less than 8" . The flux at 850 was 
determined as the integrated flux over the source extension. 



* We follow the definition of IR luminosity of Sanders & Mirabel 
( 1996 ) with Lir = L(8 - 100G>m) 
5.I6-F25 + 2.58-F60 + i^ioo) where the fluxes are in Jy. 



5.6 10 5 £^(13.48^12 
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Table 1. FIR and submillimctre data 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Name 


Fl2/jm 


F25nm 


F60Mm 


FlOO/jm 


F 8 50 M m (a) 




[Jy] 


[Jy] 


[Jy] 


[Jy] 


[mjy] 



NGC 1614 
NGC 3110 
NGC 4194 
NGC 4418 
NGC 5135 
NGC 5256- 
NGC 1653 
NGC 5936 
NGC 6240 
Arp 193 
Arp 220 
Mrk 231 
Mrk 273 
Zw 049 
IC 2810 
NGC 1667 
NGC 2623 
IRAS 0519 
IRAS 1212 



1.44 
0.59 
0.83 
0.93 
0.64 
0,23 



7.29 
1.04 
4.31 
9.32 
2.40 
0,98 



32.3 
10.7 
21.4 
40.7 
16.9 
7,3 



32.7 
19.2 
25.9 
32.8 
28.6 
11.1 



110 ± 21< 3 
144 ± 27< 2 
113 ± 22^ 
290 ± 13( x 
182 ± 20^ 
71 ± 7' 1 ) 



0.70 



1.29 



11.0 



_2H3 176 ± 24< : 



0.48 
0.56 
0.26 
0.48 
1.87 
0.24 
0.08 
0.19 
0.43 
0.21 
0.73 
0.11 



1.26 
3.42 
1.36 
7.90 
8.66 
2.28 
0.77 
0.58 
0.68 
1.74 
3.44 
0.52 



8.5 
22.7 
15.4 
112.0 
32.0 
21.7 
20.8 

5.9 

5.9 
23.1 
13.7 

8.5 



16.1 
27.8 
25.2 
114.0 
30.3 
21.4 
29.4 
10.3 
14.7 
27.9 
11.4 
10.0 



147 ± 32^ 
168 ± lO' 3 
104 ± 16' 1 
792 ± 26' 1 
126 ± 13^ 
104 ± lO' 1 
180 ± 16C 1 

< 52( 2 ) 

< 66< 3 ) 
<198< 3 ) 

<ll (i. 3 ) 

< 52« 



arguments suggest values of (3 between 1 and 2 (see Tielens 
& Allamandola 1987). Models of interstellar dust that are 
both astronomically realistic in composition, and able to re- 
produce the observed dust extinction and emission curves, 
predict /3 ~ 2 (e.g. Draine & Lee 1984, Ossenkopf & Henning 
1994, Kriigel & Siebenmorgen 1994) . The value of /3 becomes 
somewhat lower (j3 ~ 1.5, Desert, Boulanger & Puget 1990) 
if stochastic heating of small grains to high temperatures 
is taken into account. Even smaller values of /3 ~ 0.4 are 
predicted/observed in the atmospheres of Vega-type stars, 
where large dust grains that do not satisfy the cond ition of 
A>a even at submill imetre wavelengths can form (Kriigel 



Siebenmorgen 1994) 



( a ) The quoted flux error is estimated from the noise across the 

SCUBA iiia,ps, <Tn- The upper limits a,rc 3 (Th- 

The chop throws used are (1) 60", (2) 100" and (3) 120". 

4 DUST EMISSION 

We have fitted simple models to the observed 850 /im flux 
densities together with the IRAS flux densities at 100 /im, 
60 /im, and, in the case of the two-temperature model, addi- 
tionally at 25 /im. The fits are made with only three or four 
data points and so it is not possible to constrain a model 
with more than this number of free parameters. Therefore, 
a number of assumptions are made about the dust. In the 
first model we assume that the dust is at a single physical 
temperature, and that the emission is optically thin. In the 
second model we fixed the dust emissivity law, and adopted a 
two-component optically thin dust model, while in the third, 
we explored an optically thick model. In each case the data 
fits were made by varying the free model parameters and 
searching for the minimized \ 2 - 



4.1 Optically thin, single-temperature emission 

We have used a single-temperature, optically thin dust 
model given by 



Sx = Bx(T)t x oc B x (T)k x 



(1) 



to fit the 850, 100 and 60 fj,m flux densities, where Bx (T) is 
the Planck function, tx is the dust opacity and kx the dust 
absorption coefficient. At wavelengths long compared to the 
physical grain size (A 2> a, with a the grain size), the dust 
absorption coefficient can be parameterized by 

-8 



kx oc X 



(2) 



The wavelength dependence of kx is not well known, and de- 
pends on the physical composition of the dust. Theoretical 



In Fig. 1 we show the data for the galaxies together 
with the best model fits. Three different curves are shown: 
the solid line depicts the dust emission model with the best- 
fit /3 and T while the dotted and dashed lines show the 
best-fit T models with /3 of 1 and 2 respectively. The dust- 
temperatures derived for the case of variable f3 lie between 30 
and 60 K. When only upper limits at 850 ^m were measured, 
the fits were made using a fixed /3 and the 60 and 100 fim flux 
densities only. The results for (5 and T are not significantly 
different if an alternate wavelength dependence of kx is used 
at short wavelengths (kx oc A -1 for A < 250/im, Hildebrand 
1983). 

The dust emission from all galaxies is well-fit by the 
single-temperature dust model, with 1 ^ /3 < 2. The distri- 
bution of the values of the best-fit f3 as a function of Lib./ Lb 
are shown in Fig. 2, where the error bars on /? have been es- 
timated by performing fits with the 850 /im flux plus/minus 
a. The ratio Lib/ Lb is frequently used as an indicator for 
the star formation activity since Lb and Lib reflect the av- 
erage star formation rates during the last ~ 3 10 9 and ~ 10 s 
yr, respectively. Starburst galaxies have a much higher ratio 
(Lib./ Lb — 1 — 100) than normal galaxies (Lib/Lb < I). The 
interpretation of Lir /Lb as an indicator of star formation 
has to be made with caution, as the ratio is also affected by 
the amount of dust present (Calzetti et al. 1995). No trend 
of /3 with Lib/ Lb is seen. The majority of galaxies in the 
sample have values of /3 around 1.5 or higher, in agreement 
with the predictions of astronomical models (see above). The 
only exceptions are of Arp 220 and NGC 4418 which have a 
lower j3, close to 1. 

These low values of /? could be an indication of very dif- 
ferent dust properties, for example, a higher proportion of 
very small or very large grains. Although we cannot exclude 
this possibility, it is not clear why these two galaxies, that 
at face value do not differ substantially from the other mem- 
bers of the sample, should possess dust grains with different 
properties. A more plausible explanation is that some of the 
simple assumptions on which our fits are based are not cor- 
rect. In the following subsection we explore the validity of 
these assumptions. 



4.2 Optically thin, multi-temperature emission 

The modeling of the previous section was based on the as- 
sumption that all the dust was at a single dust temperature. 
This is quite simplistic as in principle the dust will be at a 
range of temperatures that depends on the local radiation 
field. Observations of spiral galaxies have shown that indeed 
at least two components, a warm and a cold one (T 20K), 
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10 1 10 2 10 3 10 1 10 2 10 3 10 1 10 2 10 3 

A [yi/.m] A [jiim] A [/xm] 



Figure 1. The dust emission spectrum for the galaxies. The filled squares are our 850 fim measurements, the crosses are IRAS data 
and the open triangle data taken from the literature (Braine & Dumke 1998; Carico ct al. 1992; Downes & Solomon 1998; Klaas et al. 
1997; Lisenfeld ct al. 1996; Rigopoulou et al. 1996; Roche & Chandler 1993; Scoville et al. 1991). The fits are done taking into account 
the data at 60, 100 and, in case of detection, 850 /im. The dashed line is a fit with (3 = 2, the dotted line for (3 = 1 and in the full line 
the fit is performed also for /3. The best-fit value for j3 and for the dust temperature T in this case are written in the graphs. 
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Figure 1 — continued 




Figure 2. The exponent of the wavelength dependence of the 
extinction coefficient, j3, is shown as a function of Ljj^/Lb- 



reasonably well-fit for these as well as other galaxies not 
shown. The discrepancy between the model and some of the 
data poi nts, e.g. with the ISO fluxes at 150 and 200 /im of 



Arp 220 (Klaas et al. 1997), could easily be improved using 
more temperature components. 

The high 850/60 and 850/100 fim flux ratios observed 
in Arp 220 and NGC 4418 necessitate the presence of large 
amounts of cold dust. The temperature of the second dust 
component needed to model the emission from these two 
galaxies is T c ~ 21 — 22, rather colder than that typically 
found for active galaxies (T = 33 K, Chini et al. 1995). In 
a more realistic multi-temperature model using more than 
2 temperature components, the temperature of the coldest 
dust component would be even lower. This is surprising as, 
a priori, it is unlikely that Arp 220 and NGC 4418 have 
colder dust than other galaxies in the sample, particularly 
given that they show among the highest values of -Lir/Lco^ 



4.3 Optically thick dust emission 

As discussed in the introduction, there is mounting evidence 
to suggest that the dust opacity in LIRGs is high. To inves- 
tigate the effects of dust opacity, the data at 60, 100 and 



are necessary to explain the dust emission spectrum from 
FIR to millimetre wavelengths (e.g. Guelin et al. 1993, Chini 
et al. 1995, Alton et al. 1998). As a first approximation to a 
more realistic model, we have applied a two-temperature fit 
to the data, with (3 — 2. 

In Fig. 3 we show the results of a two-temperature 
model fit to the 25, 60, 100 and 850 fim fluxes for a selection 
of the galaxies in the sample: Arp 220 and NGC 4418 have 
been taken as examples of galaxies with low values of best-fit 
f3, while NGC 6240 and Arp 193 are examples with higher 
values of (3 (1.4 and 1.9, respectively). The data points are 



§ The ratio Lir/Lco can De interpreted as a rough measure 
of the dust temperature using the following argument: Lqo is 
proportional to the brightness tempertature, Tf,, of the molec- 
ular gas in the galaxy. If we assume that Tf, = T, the tem- 
perature of the dust, and take Lqq as a measure of both the 
gas and dust mass (assuming a fixed gas-to-dust ratio), we ob- 
tain after reordering equation 7 and integrating over wavelengths, 
Lm/Lco <x J B\(T)k\d\/T oc T 5 , for optically thin dust emis- 
sion with (3 = 2. In contrast, if the dust is optically thick, then 
this relation changes to Lir/Lco oc J B\(T)d\/T oc T 3 . In both 
cases, a high Lm/Lco is indicative of a high dust temperature. 
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10 1 10 2 1 3 10 1 10 2 10 3 

A [/um] A [/im] 



Figure 3. The dust emission fitted with a two-temperature model for j3 = 2. The temperature of the cold (T c ) and of the warm (T w ) 
dust component are indicated. 



850 /im were fit using a dust emission model with a finite 
optical depth given by: 



Sx = B X (T) (l - e"^) , 
where 

t\ = riooMm(100/im/A)' 3 . 



(3) 



(4) 



The model fits for the two galaxies Arp 220 and NGC 4418 
are shown in Fig. 4, where the dust temperature and opacity 
have been optimized for fixed values of /3 (2, 1.8, 1.6). For 
(3 — 2, the derived dust opacities are high and in Arp 220 
the fit is only marginally consistent with the observed 150 
and 200 fim fluxes. The fit improves with decreasing (3, with 
the derived opacities becoming consistent with independent 
estimates of other authors. 

We conclude that the high 850/60 and 850/100 [im flux 
ratios of Arp 220 and NGC 4418 relative to the rest of the 
sample is consistent with a high dust opacity. 

The models discussed are very simplistic, and the true 
physical situation is most likely a combination of both a 
multiple-component dust temperature and a high dust opac- 
ity: a high dust opacity would result in the dust being 
shielded from the heating radiation, resulting in cooler dust 
components in regions further from the radiation source. 



5 DUST AND GAS MASSES 

The second objective of the project was to compare the 
molecular gas masses of the sample derived using the two 
standard methods (a) from CO measurements (b) from dust 
measurements. 



5.1 Molecular gas mass derived from the CO 
emission 

CO luminosities were calculated using CO (1-0) measure- 
ments taken from the literature. From these, we derived 
molecular gas masses using: 

Mco =4.78Lco[M ] 



(5) 

( [Sanders et al. 1991 ) , with Leo , the CO luminosity in K km 
s~ pc , given by 



Leo = 26.557} 2 /coei[Kkms~ 1 pc 



(6) 



where 7co is the integrated main beam temperature in K, D 
is the source distance in Mpc, and 0| the half power beam 
width (HPBW) of the beam in arcsec 2 . The quantity 7co 
has been taken from Sanders et al. ( 199l[ )P1 for the major- 
ity of galaxies, however, in the case of Arp 220, Arp 193, 
Mrk 2 31, Mrk 273 and NGC 6240 data from Solomon et al. 
(1997) (taken with IRAM 30m telescope) were used instead. 
The fluxes measured by Solomon et al. (1997) for three of 
these galaxies (Arp 220, Mrk 231 and Arp 193) are a facto r 
of two higher than those measured by Sanders et. al (1991). 
The origin of this difference is unclea r. We have adopted 
the values measured by Solomon et al. (1997) because of the 
higher signal-to-noise and better overall quality of their ob- 
served spectra. The values for Leo and Mco are tabulated 
in columns 4 and 7 of Table 2 respectively. The error in Mco 
is determined by the error in 7co, which we have estimated 



™ We have transformed the temperature scale used in 
Sand ers et al. (4991*1 to T Lu by multiplication with r; moon /r/B = 
1.19 (Sanders et al. 1991) 
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Figure 4. The dust emission, taking into account a finite dust opacity, according to eqs. 3 and 4, is shown for Arp 220 and NGC 6240 for 
different values of (3. The best-fit value for the dust opacity at 100 ^m, rioo^m, and the corresponding dust temperature T are indicated. 



by co mparing the data used here with other published CO 
data (Young et al. 1995), where we find CO(l-O) fluxes for 
11 galaxies of our sample. On average, there is a difference 
of about 25 per cent between the two data sets which we 
adopt as the uncertainty in 7co- 



5.2 Molecular gas mass derived from the dust 
emission 

The molecular gas mass can be calculated from the 850 [im 
emission: 



M850/jm = 



D 2 Sg5Qfj,m 



n 



(7) 



where optically thin dust emission at 850 [im is assumed 
(this is correct even for the high dust opacities estimated in 
the previous section). 

We have adopted temperatures derived in the best-fit j3 
models. We have used the value for fesso/jm given by Hilde- 
brand (1983) which, for (3 — 2 and for a gas-to-dust mass 
ratio of 150, is fegsojtm = 0.0058 cm 2 g _1 . For comparison 
with other authors, this value can be extrapolated to 1.3 
mm, where we obtain (for /J = 2) fci.3 mm = 0.0025 cm 2 g _1 . 
Based on extensive simulations of dust emission and ab- 



sorption Kriigel, Steppe & Chini ( 199C) est imate fci.3 mm = 
0.003 cm 2 g _1 . Mezger, Wink & Zylka ( |1990j ) suggest a value 
of fci.3mm = 0.0013 b z cm 2 g _1 , where b is a measure of the 
dust environment and z is the metallicity. Adopting a value 
of b — 1.9 for moderately dense gas yields (for z = 1) a 
value of fci.3 mm = 0.0025 cm 2 g _1 , which matches the value 
adopted in this paper. In contrast, the "standard" dust 



model of Draine & Lee ( 1984 ) gives a somewhat lower value 



of fci.3 mm = 0.0015 cm 2 g _1 , which has been found to agree 
with observational comparisons between dust extinction at 
IR wavelength and dust emission at FIR/mm wavelength in 
dense Galactic molecular cloud cores (Kramer et al. 1998, 
Lehtinen et al. 1998). 

The most significant uncertainty in this method, apart 
from kg^nm, is the dust temperature. The error in the dust 
temperature can be estimated by comparing the tempera- 
tures derived in the different models (optically thin/thick, 
single temperature or various temperature components) 
used here, which differ by about 50 per cent. At submillime- 
tre wavelengths the Planck function depends linearly on the 
dust temperature, so this uncertainty in the dust tempera- 
ture produces an error in Msso^m of 50 per cent also. Adding 
to this quadratically the measurement error of the 850 ^m 
flux densities (on average 25 per cent), we estimate an error 
of 55 per cent for Mgso^m- 



5.3 Comparison between the two masses 

In Fig. 5 we show the comparison between the molecular gas 
mass derived from the CO emission, Mco, and the molecu- 
lar gas mass derived from the 850 /im continuum, Afeo^m- 
Mco is systematically higher than Msso^m ■ A good correla- 
tion (correlation coefficient r — 0.82) exists between the two 
quantities, spanning nearly 2 orders of magnitude. 

In Fig. 6 we have plotted the ratio between the molecu- 
lar gas mass derived using the two alternative methods as a 
function of Lir and as a function of Lir/Lb to see whether 
there is any trend with luminosity or the star formation ac- 
tivity. The error of Mco/Aisso^m is the root of the sum of 
the squares of the errors in Mco and M85o M m, and so ap- 
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Table 2. Luminosities and gas masses 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


Name 


Distance 


Lm, 


Leo 


Lb 




M C o 




L B 


J'co 
8 M m 




[Mpc] 


[lO n L ] 


[10 s L' ] 


[1O 1O L ] 


[1O 9 M ] 


[10 9 M Q ] 








NGC 1614 


62.3 


3.8 


31.9 


1.7 


2.9 


15.3 


117.8 


22.2 


5.2 


NGC 3110 


69.3 


1.6 


59.7 


1.7 


5.9 


28.6 


27.1 


9.4 


4.8 


NGC 4194 


41.4 


1.1 


6.7 


1.1 


1.4 


3.2 


163.6 


10.0 


2.3 


NGC 4418 


34.2 


1.3 


5.3 


0.4 


1.6 


2.5 


247.2 


34.4 


1.6 


NGC 5135 


56.9 


1.7 


34.4 


3.3 


5.2 


16.4 


49.2 


5.2 


3.2 


NGC 5256 


115.9 


2.9 


63.0 


3.3 


7.6 


30.1 


45.6 


8.7 


3.9 


NGC 5653 


54.4 


1.1 


15.0 


2.3 


4.5 


7.2 


72.1 


4.7 


1.6 


NGC 5936 


59.0 


1.0 


19.0 


2.2 


4.4 


9.1 


52.4 


4.5 


2.0 


NGC 6240 


100.9 


6.1 


90.4 


4.9 


12.1 


43.2 


70.4 


12.9 


3.6 


Arp 193 


97.8 


4.0 


44.3 


1.6 


9.1 


21.2 


91.3 


25.2 


2.3 


Arp 220 


79.2 


15.8 


88.0 


2.0 


30.0 


42.0 


179.9 


78.2 


1.4 


Mrk 231 


173.9 


31.0 


85.6 


7.0 


24.2 


40.9 


361.8 


44.2 


1.7 


Mrk 273 


153.8 


12.2 


57.8 


5.1 


15.3 


27.6 


211.8 


24.1 


1.8 


Zw 049 


52.3 


1.3 


8.9 


0.2 


3.8 


4.3 


149.9 


81.5 


1.1 


IC 2810 


140.5 


3.2 


54.8 


3.4 


< 10.0 


26.2 


57.6 


9.3 


> 2.6 


NGC 1667 


58.9 


0.8 


12.3 


3.4 


< 2.5 


5.9 


62.3 


2.3 


> 2.3 


NGC 2623 


76.1 


3.2 


26.6 


2.8 


< 9.4 


12.7 


121.1 


11.7 


> 1.4 


IRAS 0519 


167.5 


11.7 


69.8 




< 21.2 


33.4 


167.4 




> 1.6 


IRAS 1212 


293.3 


16.7 


123.5 




< 36.7 


59.0 


135.1 




> 1.6 



,, : 75 km s 1 Mpc 1 and taking into accoun t the Virgocentric flow. 

(3) L IR = L(8 - lOOO^m) = 5.6 10 5 1)2^(13.48^12 + 5.16F 25 + 2.58F 60 + Fioo) where the flux is in Jy ( ganders fc Mirabel 1996] ). 



(2) Distance taken from Sanders et al. (1991), based on Hq 

'Mpc 



(4) L'=K km s" 

(5) Calculate d from B^,, take n from de Vaucouleur et al. (1991), except for Zw 049 and NGC 3110 which have been taken from Lehnert 
& Heckman (199J) and Lutz ( 1992j ), respectively. 

(6) Molecular gas mass estimated from the CO emission, calculated according to eq. (5). 

(7) Molecular gas mass estimated from the dust emission, calculated according to eq. (7). 



o 



ID 

CO 




M C o [M ] 

Figure 5. The molecular gas mass derived from the CO- 
luminosity, Mqq, and the molecular gas mass derived from the 
dust emission, Mg5o Mm . The diagonal line shows the relation 
Mqq = Mg50 Mm . 



proximately 60 per cent. The average ratio, Mco/Mg 5 o^m is 
2.6 ± 1.3, ignoring upper limits. There is no trend visible for 
Mco /Msso/im to vary with Lm- There might be a weak anti- 
correlation in the LIRG sample between Mco/A^850/jm and 



Lir/Lb, but the statistical significance is low and a larger 
sample size will be needed in order to confirm or refute this 
anticorrelation. 

Also plotted in Fig. 6 are data taken from the litera- 
ture for other galaxies, including spiral galaxies (NGC 3637, 
NGC 891) and actively star forming and starburst galaxies 
(Markarian galaxies, NGC 253). The Mco/A/ssoum ratios 
for these comparison galaxies have been rederived using the 
same parameters as for the LIRG sample. The most signif- 
icant differences between the LIRG sample and these other 
galaxies is seen in the spiral galaxies NGC 891 (factor of 3.3) 
The starburst galaxy NGC 253 has about the same ratio as 
the average of the LIRGs, while the Markarian galaxies have 
a lower value, possibly due in part to the fact that Mqq has 
been derived from the CO(2-l) line emission assuming a ra- 
tio of CO(2-1)/CO(1-0) of 1. 

We have found for the LIRG sample that the molec- 
ular gas mass derived from the CO is higher than that de- 
rived from the 850 fitn. This confirms the conclusion of other 
authors, that the standard galactic CO-to-molecular mass 
conversion factor overestimates the molecular gas mass in 
LIRGs. The resultant ratio of Mco/-^f850/jm is slightly less 
than that seen in other studies, the results of which sug- 
gest that the CO emission overestimates the molecular gas 
mass by a factor of about 4-5 (Solomon et al. 1997, Downes 
& Solomon 1998) or more (factor 3-10, Shier et al. 1994, 
Braine & Dumke 1998). If this higher value is true, it im- 
plies that the molecular gas mass derived from the 850 fim 
flux is also an overestimate. This could be possibly due to a 
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Figure 6. The ratio between the molecular gas mass derived from the CO emission and the molecular gas mass derived from the dust 
emission as a function of Ljjj (a) and of Ltr/Lb (b) for the LIRG sample (filled squares and arrows). For comparison, data from the 
literature (crosses) for spiral galaxies (NGC 3627, Sievers et al. 1994; NGC 891, Guelin et al. 1993), Markarian galaxies ("Mrk", the error 
bar gives the range for the values in the sample; Chini, Kriigel & Steppe 1992), and the starburst galaxy NGC 253 (the error bars give 
the range of values found in the center and in the disk; Mauersberger et al. 1996) are also included. 



higher dust-to-gas mass ratio in these galaxies which could 
be caused by the higher metallicity produced by intense star 
formation as seen in galac tic nuclei and bulg es, for example 
in the inner p art of M 51 (|Guelin et al. 1995 ) and the centre 
of NGC 253 flMauersberger et al. 1996^ 
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6 CONCLUSIONS AND SUMMARY 

We have presented new data at 850 /im taken with SCUBA 
at the JCMT for 19 luminous infrared galaxies. Fourteen 
galaxies were detected at a level of 5cr or better. Using the 
850^im data as well as IRAS fluxes taken from the literature, 
we have shown that the submillimetre/far- infrared spectral 
energy distribution for the majority of the sources in the 
sample is well-described by an optically thin, single temper- 
ature dust model, with /3 ~ 1.4 — 2. A lower value of j3 ~ 1 
is required to fit the data for two of the galaxies, Arp 220 
and NGC 4418. 

This difference can be attributed to one, or a combina- 
tion of: (i) different intrinsic dust properties, (ii) the pres- 
ence of a large cold-dust component, (iii) high dust opacities. 
We have discussed the various possibilities and conclude that 
the most physical explanation is a high dust opacity. 

We have compared the molecular gas masses derived 
from the dust emission to the molecular gas masses derived 
from the CO emission of our sample. We find, on average, 
that Mco is a factor of 2-3 higher than Mg5o Mm . If the CO 
luminosity indeed overestimates the gas mass by a factor 
of about 5-10 as indicated by other studies, then the dust 
emission must also overestimate the gas mass. The most 
likely reason for this is a higher dust-to-gas mass ratio in 
LIRGs. 
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